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Abstract This tutorial article gives an introduction to the methods needed to 
treat interacting electrons in a quantum wire with a single occupied 
band. Since one-dimensional Fermions cannot be described in terms 
of noninteracting quasiparticles, the Tomonaga-Luttinger model is pre- 
sented in some detail with an emphasis on transport properties. To 
achieve a self-contained presentation, the Bosonization technique for 
one-dimensional Fermions is developed, accentuating features relevant 
for nonequilibrium systems. The screening of an impurity in the wire is 
discussed, and the insight gained on the electrostatics of a quantum wire 
is used to describe the coupling to Fermi-liquid reservoirs. These parts 
of the article should be readily accessible to students with a background 
in quantum mechanics including second quantization. To illustrate the 
usefulness of the methods presented, the current-voltage relation is de- 
termined exactly for a spin-polarized quantum wire with a particular 
value of the interaction parameter. This part requires familiarity with 
path integral techniques and connects with the current literature. 

Keywords: Tomonaga-Luttinger liquid, Bosonization, Electronic transport proper- 
ties in one dimension, impurity scattering, current-voltage relation. 
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1. INTRODUCTION 

Within the last decade there has been increased interest in the be- 
havior of quasi one-dimensional Fermionic systems, due to significant 
advances in the fabrication of single channel quantum wires [1, 2, 3] 
based on semiconductor heterostructures and the observation of non- 
Fermi liquid behavior in carbon nanotubes [4, 5, 6]. While the unusual 
equilibrium properties of Fermions in one dimension have been studied 
since many decades and are well documented in review articles [7, 8], 
nonequilibrium quantum wires are an area of active research with many 
important question remaining to be answered. 

In this article, we give a rather elementary introduction to the theoret- 
ical framework underlying much of the present studies on transport prop- 
erties of one-dimensional Fermions. While we do not review extensively 
features of the Tomonaga-Luttinger model [9, 10] upon which these stud- 
ies are based, we give an elementary introduction to the Bosonization 
technique which is an essential ingredient of current theoretical methods. 
We do not review the rather long history of Bosonization starting with 
the work by Schotte and Schotte [11] in 1969. Some important articles 
are contained in a book of reprints collected by Stone [12]. Our approach 
is based on Haldane's algebraic Bosonization [13], which can be under- 
stood with the usual graduate level background in physics. For a more 
in-depth discussion of the method, we refer to a recent review by von 
Delft and Schoeller [14]. The field theoretical approach to Bosonization, 
which is probably harder to learn but easier to apply, has lately been 
expounded by Gogolin, Nersesyan and Tsvelik [15]. 

We employ the Bosonization technique to describe a quantum wire 
coupled to Fermi-liquid reservoirs. In this connection the electrostatic 
properties of the wire play an important role. Landauer's approach [16] 
to transport in mesoscopic systems, which is based on Fermi liquid the- 
ory, is generalized to take the electronic correlations in a single channel 
quantum wire into account. In the Bosonized version of the model the 
coupling to reservoirs is shown to be described in terms of radiative 
boundary conditions [17]. This allows us to use the powerful Bosoniza- 
tion method also for nonequilibrium wires. About a decade ago, Kane 
and Fisher [18] have noted that transport properties of one-dimensional 
Fermions are strongly affected by impurities. Even weak impurities have 
a dramatic effect at sufficiently low temperatures leading to a zero bias 
anomaly of the conductance. It is the aim of the article to present the 
theoretical background necessary to study the recent literature on this 
subject. Again, we do not provide a review of transport properties of 
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the Tomonaga-Luttinger model. Rather, the methods developed are 
illustrated by treating a particular case. 

1.1. NONINTERACTING ELECTRONS IN 

ONE DIMENSION 

Let us start by considering first noninteracting electrons of mass m 
moving along a one-dimensional wire with a scatterer at a; = 0. For 
simplicity, the scattering potential is taken as a J-potential 

V,,{x) = -K5{x) (1) 

m 

where A characterizes the strength. The Schrodinger equation 

i)"{x) + —K5ix)Mx)=ei){x) (2) 

2m m 



has for all positive energies 



a solution (fc > 0) 



efc = 1^ — (3) 

Zm 



describing a wave incident from the left that is partially transmitted and 
partially reflected. The transmission amplitude 

determines the transmission coefficient T^ and the phase shift rj^. Like- 
wise, there is a solution {k > 0) 

V'-fc(x) = ^k{-x) (6) 

describing a wave incident from the right. 

When the ends of the wire are connected to electrodes with electrical 
potentials ^l and /i/j, a voltage 

U = ifiL- f^R)/e (7) 

is applied to the wire, where e is the electron charge, and an electrical 
current 

/ + 00 
dkfkjk (8) 

-oo 



flows. Here 

jk (x) = Im V-fc (x) —-K-i^k (x) = Tk (9) 

m ox lirm 

is the particle current in state ipk which is independent of x, and the 

f _ / fi^k- f^h) ,k>0 , . 

with /(e) = \/{e^^ + 1) are state occupation probabilities determined 
by the Fermi function of the electrode from which the particles come. 
Both electrodes are assumed to be at the same inverse temperature f3. 
Putting ixl,r = ef i \eU ^ we readily find for small voltages U 



GU ill) 



with the conductance 



G = e' dk^n^fisk-SF) 
Jo zirm oep 



(12) 



p2 (-oo 



-|/(£-eF! 



Provided (3eF is large, this yields the Landauer formula for a single 
transport channel 



e2 



G=jTf, (13) 

where Tp is the transmission coefficient at the Fermi energy ef- If we 
take into account the spin degeneracy of real electrons, the conductance 
becomes multiplied by 2. 

1.2. FANO-ANDERSON MODEL 

Of course, to describe electrons in one dimension, we may also start 
from a tight binding model with localized electronic states at positions 
Xj = aj where a is the lattice constant. The Hamiltonian in the presence 
of an impurity at x = then takes the form of the Fano- Anderson model 
[19, 20] 

+ 00 

j=-oo 

where t is the hopping matrix element and Sq > is the extra energy 
needed to occupy the perturbed site at x = 0. The aj are Fermi oper- 
ators obeying the usual anti-commutation relations. This model is also 
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exactly solvable [21] with eigenstates in the energy band 

Ek = —2tcos{ka) ,\k\ < — , (15) 

a 

and the transmission amplitude takes the form 

1 

1 + i/\ 



tk = rr—TV- ' (16) 



where 

Xk = — sin(A;a) = . (17) 

The last relation follows by virtue of v^ = {l/fi)dek/dk. It is now 
easily checked that this model leads to the same result for the linear 
conductance than the free electron model examined previously, provided 
we match parameters of the unperturbed models such that the Fermi 
velocities vp at the two Fermi points coincide, and we adjust the strength 
of the (5-function in Eq. (1) such that 

h^ 

— A = eoa. (18) 

m 

Then, we end up with the same transition coefficient Tp at the Fermi 
energy. 

1.3. NONINTERACTING 

TOMONAGA-LUTTINGER MODEL 

Apparently, for low temperatures and small applied voltages, the con- 
ductance only depends on properties of states in the vicinity of the Fermi 
energy ep. We can take advantage of this fact by introducing still an- 
other model, the noninteracting Tomonaga-Luttinger (TL) model, which 
has the same properties near the two Fermi points zizkp but is more 
convenient once we introduce electron interactions. Let us formally de- 
compose the true energy dispersion curve Sk into two branches + and 
— of right-moving and left-moving electrons, respectively, where these 
branches comprise states in the energy interval [ep — A,ep + A], (cf. 
Fig. 1.1). If A is chosen large enough, these branches should suffice 
to describe the low energy physics of the true physical model, since for 
low temperatures and small applied voltages states with energy below 
ep — A are always occupied while states above ep + A are empty. The 
Hamiltonian of the noninteracting TL model reads 

H = J2J2 ^P.fc [4,fcS,fc - i4>,k^P,k)o] ■ (19) 

P k 



\ 


ep+A 


^ 


\ 


k '^ 




h 




/ 




\^-A 




> 



Figure 1.1. Physical energy dispersion curve ek (thin hne) and the two branches 
+ , — (thick lines) of the TL model 



Here j? = it labels the two branches. We have introduced a quantization 
length L such that wave vectors are discrete^ 



k 



27r 

T 



nfc , Hk integer. 



(20) 



Further, the £±{k) are single particle energies measured relative to the 
Fermi energy ep, i.e., £j^{kp) = e-{—kF) = 0. The operators c , 
and Cp^k are Fermi creation and annihilation operators obeying anti- 
commutation relations, in particular 



'^p,k ' '^p'k 



Sp,p'h,k' 



(21) 



The sum over k states in Eq. (19) is restricted to k values near zizkp such 
that e±{k) £ [—A , A]. Finally, the ground state energy is subtracted in 
Eq. (19) where the ground state |0,0)o is defined by 



c+ fc|0, 0)o = c_ _fc|0, 0)o = ioT k> kp , 



(22) 



c^+,fc|0,0)o = cl _JO,0)o = iorkKkp. 
Frequently, the spectra e±^k are linearized about the Fermi points, i.e., 

e±^k = ±fiVF{k =F kp) , (23) 



^We use periodic boundary conditions here and consider the limit of large L in the sequel. 
The same techniques can also be used for other boundary conditions [14, 22]. 
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and then the cutoff energy A is used as a large energy scale regularizing 
divergent expressions. Since for low energy systems only inert empty 
or occupied states are added, an increase of A is admissible. However, 
the linearization of the spectra is only realistic in the close vicinity of 
the Fermi points. In fact, some phenomena not discussed here, e.g. the 
thermopower [23], depend on band curvature. We shall come back to 
the limitations of the linearization (23) below. 

Before we discuss transport properties of the TL model, we first in- 
troduce methods, the advantage of which becomes apparent only when 
we pass on to the case of interacting electrons. These methods are inde- 
pendent of the precise dispersion law as long as the Ep^k ^-re monotonous 
functions of k. 

2. BOSONIZATION 

The noninteracting TL model allows for a formulation in terms of 
Bose operators. We discuss the Bosonization technique here for spinless 
Fermions first, and then extend it to the spinful case. 

2.1. DENSITY OPERATORS AND THEIR 

ALGEBRA 

While the physical problem of one-dimensional Fermions is described 
by a single energy dispersion curve £k with empty states at both ends of 
the range of k values, the TL model introduces two branches with empty 
states at one end but occupied states at the other end of the k range. 
This leads to unusual algebraic properties we will discuss now. Since the 
range of allowed k values plays an important role in this discussion, we 
keep track of it in detail by introducing 

Let us define Fourier components of the densities of p-movers {p = it) 

by 

Pp,q = Z] Wp,kWp,k+q cl^kCp,k+q ■ (25) 

k 

Since Fermi operators on diff'erent branches anti-commute, we have 

[Pp,q , Pp',q'] - = for p / p . (26) 

On the other hand, for the commutator on the same, say the + branch, 
we find 

[p+,q,P+,q']- = Y. ^+<kW+,k+qW+^k'W+^k'+q' (27) 

kk' 



'+,k^+,k+q(^+,k'^k'+q 



'^+,k'^+,k'+q'(^+,k'^+,k+q 



We now use anti-commutation relations for the second and third opera- 
tor in each of the two products of four Fermi operators, e.g., c , 



,k+q'^+,k' 



^k+q,k' — c\- k'^+ k+q- Then, the remaining products of four Fermi opera- 
tors are easily seen to cancel by virtue of the anti-commutation relations, 
and the two terms with two Fermi operators can be written as 



[P+,q-, P+a']- 
k 



(28) 



Wl 



,k+q 



Wl 



,k+q' 



H,k^+Mg+(]' ■ 



The terms of the sum are finite only for W^^k = M^+.fc+g+g' = 1- But 
then W^^k+q and W^^k+q' are also equal to 1 if (7 and q' have the same 
sign. In this case all terms of the sum (28) vanish. Hence, a nontrivial 
commutator can only arise if q and q' have different signs, say q > and 
q' < 0. Then non-vanishing terms of the sum (28) may occur near the 
upper cutoff where e+ ^ ~ A, in particular, when 



W. 



W^ 



+ ,k+q+q' 



w_ 



+ ,k+q' 



1 , but W^ 



+ ,k+q 



0, 



as illustrated in Fig. 1.2. However, in this case the operator c 



+,k'-+,k+q+q' 

tries to annihilate a particle in the state k + q + q' with an energy near 
A. In the low energy sector of the model these states are always empty 
and c^^k+q+q' can be replaced by zero. Hence, for a low energy system 
we get no contribution to the commutator from states near the upper 
cutoff energy. 

A contribution near the lower cutoff energy may arise from terms with 
W+^k = W+^k+q+q' = W+^k+q = 1 but W+^k+q' = siuce e+,fc+g' < "A. 
In this case we use the anti-commutation relations to write Eq. (28) as 



[p^ 



,q ' 



P4 



k 



W+^kW+^k+q+q' 



WIm 



wi 



,k+q' 



-q' ^+,k+q+q'^+,k 



(29) 



J 



With a similar argument as above, the operator c, ^ may now be re- 
placed by zero since near the lower cutoff all states are occupied. We 
thus obtain 



[P^ 



.<? ' 



P+,q'. 



k<kp 



Wl,k 



w_ 



+ ,k+q 



Wl,k^ 



(30) 



where the condition k < kp ensures that we do not get contributions 
from the region near the upper cutoff where the Fermi operators in 
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k+q k k+q 



Figure 1.2. Wave vectors in the sum (28) for g > and g' < near the upper cutoff 
energy of the + branch. 



Eq. (29) cannot be dropped. Finally, since 



E wlkW\ 



+ ,k+q 



k<kf 



k<kp 



Wlk- 



,Wlk 



n 



Q ' 



(31) 



where q = ^Uq determines the number Uq of additional non-vanishing 
terms in the first sum, and with a similar reasoning for the — branch 
and other signs of g, we get 



[Pp,q 5 Pp',q' . 



P^p,p'5q-qinq. 



(32) 



As we have seen this nontrivial commutator only arises since we have in- 
troduced two branches with empty states at one end and occupied states 
at the other end. The result (32) can also be derived in the same way 
when the sharp cutoff functions W±^k defined in Eq. (24) are replaced by 
a smooth cutoff. Furthermore, the relation (32) is not an exact operator 
relation but holds only in the low energy sector of the Fock space. How- 
ever, for technical convenience, we may use the linearized spectrum (23) 
and send the cutoff A to infinity. We then obtain a model where the 
relations (32) hold as formally exact commutation relations but have to 
remember that only the low energy properties of this model are related 
to the physical problem of one-dimensional electrons. 
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+ branch 




Figure 1.3. 



Action of bq for q > on states of the + branch. 



2.2. BOSE OPERATORS AND BASIS 

VECTORS 

With the result (32) it is now straightforward to introduce Bose an- 
nihilation and creation operators^ 



—I 



Y^ G{pq)P: 



p,g ' 



q\ P 



(33) 



bl 



Yl ^iPl)P: 



p-q ' 



riql p 



where q j^ 0. Inserting (25), we see that for g > the operator bg 
lowers the wave vector k of right-movers by q, provided the state k — q 
is empty, cf. Fig. 1.3, while b-g acts accordingly on left-movers. By 
virtue of Eq. (32), the bg and 6j satisfy Bose commutation relations, in 
particular 



h 



i'^- 



s, 



q,q 



(34) 



When acting on a A^-particle state, the operators b' and bg create and 
annihilate electron-hole pairs, respectively, but they conserve of course 
the particle numbers Np of each branch. Consider now A^_|_ -particle 
ground states of right-movers where |A^+ = 0)o is the many-body state 
with all single particle states k < kp occupied and all states k > kp 
empty, while in the state lA*".!. = l)o also the level k = kp + ^ is 



^As usual, 9{x) denotes the step function: 9{x) = 1 for a:: > 0, 9{x) = for a:: < 0. 
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C) 
C) 
() 

<> 



C) 



empty 



occupied 



|N^=-1) |N^=0) |N^=1> 

Figure 1.4- Graphical representation of A'^+ -particle ground states of right-movers. 

occupied, and so on, as illustrated in Fig. 1.4. One can then show that 
the states 

g>0 V^l- 



liV- 



+/0 > 



(35) 



with iV+ = 0, ±1, ±2, ... and ruq = 0, 1, 2, 



form a complete basis in 
the Fock space of right -movers. It is obvious that the states (35) are all 
within the Fock space spanned by the standard basis vectors |{A^_|_^fc}), 
where the A^+,fc = 0, 1 are the usual Fermi occupation numbers, that is 
eigenvalues of c^ ^c^ ^,. Haldane [13] has demonstrated completeness of 
the basis (35) with the help of the model Hamiltonian 

1. 

2' 



h. 



E 

k 



(nk 



'^+,k'^+,k 



(36) 



where 



"+,k'^+,k 



^+,k^+,k 



for 
for 



A;>0 



(37) 



This Hamiltonian assigns positive energies to all states and allows for an 
explicit evaluation of the partition function Z = ^^ exp(— /3£'a) for the 
two sets {a} of basis vectors. Since Z is a sum of positive terms, and 
the sum over the states (35) gives the same result as the sum over the 
states |{A''+^fc}), the basis \N^,{mq}q^o) spans the entire Fock space of 
right -movers. 
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Likewise we may introduce a basis 

|iV_,{mJ,<o) = n ^^l^^-)o (38) 

in the Fock space of left -movers and combine both to the basis 

|iV+,iV_,{mJ) = n^^|iV+,iV_)o (39) 



in the Fock space of one-dimensional particles. 

2.3. LADDER AND PARTICLE NUMBER 

OPERATORS 

It is clear that the Bose operators 6„, 6j cannot generate the whole 
algebra of operators in the Fock space since they preserve the particle 
numbers Np. Hence, the Bose operators need to be supplemented by 
ladder operators removing or adding a particle of branch p. The lowering 
operators Up are defined by 

Up\Np, N.p)o = p^+ \Np - 1, N^p)o , (40) 

where the sign factor p^+ is one possible choice assuring anti-commutation 
relations for Fermi operators on different branches, and 

[Up,bg]^ = [Up,bl]^=0. (41) 

The adjoint raising operator obeys Uj, = Up^. The action of the ladder 
operator Uj^ on a state in the Fock space of right -movers is illustrated 
in Fig. 1.5. Since the ladder operators commute with the Bose operators 
6„, 6L they preserve the electron-hole pair excitations present in a state. 

Now, the particle number Np is an eigenvalue of the particle operator 
^p = H{ 4>,kCp,k - (4,fcS,fc)o) ' (42) 

k 

where the average ( )o is over the \Np = 0)o ground state. The ladder 
operators and the particle numbers satisfy the commutation relations 

Np, Up' _ = -Sp^p'Up , 

Np,ul]_ = 6p,p'UJ,, (43) 
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o empty 
• occupied 



^) 



UI^A) 



Figure 1.5. Illustration of the action of U+ on a many-body state |?/>) of right- 
movers. 

that are easily shown by applying them to an arbitrary vector from the 
basis set (39). Furthermore, the ladder operators obey 






ulul 



0, 



2 6. 



'p,p' ) 



(44) 



which may be demonstrated in the same way taking the sign factor in 
Eq. (40) into account. 



2.4. 



BOSONIC PHASE FIELDS 



With the Bose and ladder operators we should be able to represent 
all operators in the Fock space of the TL model. To demonstrate this 
one usually constructs explicitly the Fermi annihilation operators 

^p(x) = 4^5]e^'=-'cp,fc. (45) 



^/L 



It is convenient to first introduce Bosonic phase fields 



^p{x) 



1 

2^ 



E ^(P9) 



1 



g^O 



Jqx r 



9' 



and associated Hermitian phase fields 

^p{x) = ipp{x) + ifUx) 



(46) 



(47) 



14 



These are related to the densities 
Pp{x) = ^l{x)^p{x) - (^t(x)^p(x))o 



1 
I 






(48) 



where the pp^q = J2k ^p k^p k+q ^^^ ^^^ Fourier components studied above. 
Using Eq. (33), the phase field (46) can also be written as 



lfp{x) 



iE%9)u^«"^pp,.- 



27r 



<?7^0 



\nn 



(49) 



Because of the factor 9{pq), we may replace \nq\ by pUg = 2^pq- Then, 
adding the Hermitian conjugate field, we readily find for the gradient 



d 



p 






,iqx ~ 



1 



Pp.g = V Pp{x) - yNj 



L 



(50) 



The Bosonic phase fields ^{x), ^'^{x) obey the commutation relations 



[fp{x) , fp'{x')]_ = iflix) , fl,{x') 



while by virtue of Eqs. (34) and (46) the commutator 



"Ppix), fl'{x') 



Since 



E- 



^P,P' V^ }_ ip^(x-x')n 

47r2 ^ n 

n=l 



ln(l - y) , 



(51) 



(52) 



(53) 



the right hand side of Eq. (52) is logarithmically divergent for x = x' . 
However, we have to remember that the sum over n originates from a 
sum over wave vectors q = {2tt /L)n which due to the energy cutoff are 
restricted to small q. The effect of a cutoff is seen when we multiply 
the terms of the sum (52) by an exponential cutoff function exp(— ag) = 
exp(— -£^ an) which limits the sum to wave vectors q^l/a, where o is a 
cutoff length related to a cutoff energy A ~ hvp/a. We then obtain 



fpix), (Pp,[x 



'^p,p 



In 1 — exp 



27r 

T 



a — ip{x — x)] 



■47r2 



j- In f — [a — ip{x — x')'\ 



(54) 
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where the second equahty holds for L ^ a, jx — x'|. For later convenience 
we note that in particular the commutator for vanishing distance 

"Ppix) , ^l,ix^ _ = -^ In -J- (55) 

is cutoff and size dependent. 

Simpler algebraic properties are found for the Hermitian phase field 
(47). From Eqs. (51) and (54) we readily find 

[%{x) , %,{x')\ = -^ ^^ a-ip{x-x') 
I v\ j^ P\ ji- ^^2 a + ip{x-x') 

where L has dropped out, and we have for |x — x'| ^ a 

['^p{x),%,{x')]^ = !^sign(x-x'). (57) 

Here sign(x — x') is the sign of x — x' which produces a step at x = x' 
that in the presence of a cutoff is smeared over a length of order a. 
Instead of the fields $±(x) we shall mainly use the linear combinations 



(p{x) = Vvr['l'+(x) + $_(x)] 



(58) 



t9(x) = ^[^+{x)-^^{x)] 
that are readily seen to obey the commutation relations 

[<p{x) , 0(x')] _ = [t?(x) , ??(x')] _ = (59) 

and 

[(/)(x) , 'd{x')]_ = - sign(x — x') , (60) 

which shows that 

Il4x) = -h^(l){x) (61) 

is the field conjugate to t?(x) with the canonical commutator 

[i?(x) , n^(x')] = in6{x - x) . (62) 

Likewise n^(x) = —fi-^d{x) is the field canonically conjugate to (/>(x). 
Finally, we note that by virtue of Eq. (50) 

^|-^(x) =p+(x) +p_(x) - ^+ + ^- (63) 

yTT ox L 

describes the fluctuations of the total density of right- and left-movers. 
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2.5. 



BOSE REPRESENTATION OF FERMI 
OPERATORS 



To construct the Fermi operators ^p(x) in terms of the Bosonic fields, 
we start by evaluating the commutator 



[bq , ^p(a;)]_ = ap^g{x) ^p(a;) , 



where 



0{pq) e 



~iqx 



(64) 
(65) 



This follows readily by inserting Eqs. (33) and (45) into (64) with the 
help of the Fermi commutator [pp'^q-, Cp,k\_ = ~^p,p' <^p,k+q- In the same 
way we find 

bl,^p{x)]_=a;jx)^p{x). (66) 

Operating with 6q^p(x) onto a A^-particle ground state |A^-|-,A^_)o we 
have 



bg^p{x)\N+,N_)o = [bq,^p{x)]_ \N+,N_)o 

= Qp,g(x)^p(x)|iV+,iV_)o, 



(67) 



where the first relation holds because 6q|A^-(-, A^_)o = 0. Now, Eq. (67) 
shows that ^p(x)|iV+, A^_)o is an eigenstate of bg with eigenvalue ap,g(x). 
Eigenstates of Bose annihilation operators are known as coherent states, 
and from their familiar properties it follows that ^p{x)\N^,N-)o is of 
the form 



*p(x)|iV+, iV_) = Xp{x) exp Y. «m(^) H \^p - 1' ^-p)o ' (68) 

V^^o J 

where Xp{x) is an as yet undetermined normalization factor, and we have 
taken into account that ^p{x) reduces the number Np of p-movers by 1. 
With the help of the Bose commutator (34) , which implies 



bq , exp I ^ ttp^g' (x) b\, 

9Y0 



Q 



p,q{x) exp I Y^ ap^q>{x)b\, I , (69) 



it is readily seen that the ansatz (68) indeed satisfies Eq. (67). From the 
definition (46) of the Bosonic phase field fp{x) and Eq. (65) we find 



Y ap^q{x)b\ = 2TTi(fl{a 

q^O 



(70) 
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and hence from Eq. (68) 

o{Np - 1, iV_p| ^p(x) |iVp, Af_p)o = Ap(x) 
where we have made use of 



o(iVp - l,A^_p| exp 2TTi(pl{x) 



\N„ 



l,iV- 



P/O 
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(71) 



(72) 



which is easily seen by expanding the exponential. We now insert the 
representation (45) of ^p(x) into Eq. (71) and use 



N, 



)(iVp - 1, iV_p| cp,fc \Np, N^p)o = p''+ Sk,k, , 



where 



, 27r 

k„= p[ kp + —Np 



(73) 



(74) 



The sign factor p^+ arises from the anti-commutation relations of Fermi 
operators on different branches in accordance with Eq. (40). This yields 






(75) 



which can now be combined with the definition (40) of the lowering 
operator to obtain from Eq. (68) 

^p{x)\N+,N^)o = ^e^P('=f+^^^)^+27r*^^(^)f/p|iV+,iV_)o. (76) 



This represents the action of a Fermi annihilation operator on a A^- 
particle ground state in terms of Bose and ladder operators. Since 
particle-hole excitations are created by Bose operators, it is straight- 
forward to generalize Eq. (76) for arbitrary states in the Fock space. For 
a basis vector (39) we have as a consequence of the commutator (66) 



^p{x)\N+,N^,{mg}) 

\0q — Oipq\X) 



= E 

Using Eq. (76), we find 

^p(x)|iV+,iV_,{m,}) 
1 



v^ 



M'p(x)|iV+,iV_)o. 



(77) 






[bl-a;,g{3 



Vrn^- 



(78) 
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since the Bose creation operators commute with ipl{x) and Up. In view 
of 



Z^ a*p,q(.^)bg = ■ 


-2TTi 


^p{x) 










(79) 


and the famihar Bose relation e b' e~ 


■ab 


6t + a 


we 


have 








^2TTiipp{x) ^t g-27rj(/3p{x) 


= "1 


- a*p,qi 


x), 








(80) 


and Eq. (78) may be written as 
















^p{x)\N+,N^,{mg}) 














(81) 




(x) g2 


lni(pp{x) 


\N+ 


,iV_, 


{m 


J). 





Here we have made use of ipp{x)\N^, N^)o = 0, and have replaced in 
the exponent Np by the operator Np, which is appropriate since Np 
commutes with the Bose operators and thus acts on an eigenstate with 
eigenvalue Np. Now, the operator acting on the basis vector on the right 
hand side of Eq. (81) has the same form for any vector, and we thus 
obtain the operator identity 

^p{x) = ^Up e'P^^F+^ %)^ gS^i^k^) e2^*^p(^). (82) 

This is the desired representation of the Fermi operators ^p{x) in terms 
of Bose and ladder operators. The result (82) is in normal ordered form 
with all Bose creation operators to the left of the annihilation operators. 
It should be noted that this remarkable relation is independent of the 
form of the Hamiltonian, it just relates the form (45) of the Fermi oper- 
ator with an obvious meaning in the standard occupation number basis 
|{-^+,fe}) {^-,fc}) of the Fock space to a representation with an obvious 
interpretation in the basis |A^_|_, A^_, {mg}). 

For much of the following discussions another form of ^p{x) is often 
more convenient. Using the commutator (55) and the operator relation 

^A^B^^A+B^\[A,B]-^ (83) 

which holds if [^,i3]_ commutes with A and -B, we may transform the 
relation (82) to read 

^p(a;) = -^ Up e^P('=F+^ Np)x+2m'^p{x) ^ (34) 

V27ra 

where ^p{x) was introduced in Eq. (47). Finally, in terms of the phase 
fields (/>(x) and 'Q{x) defined in Eq. (58), we have 

^p(a;) = ^^ Up e»p[^>^+x NpX+V^&{x)\+iV^.p{x) ^ (35) 
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which is the form employed mostly in the hterature.^ 

2.6. BOSE REPRESENTATION OF THE 

HAMILTONIAN 

In the fohowing, we shall make use exphcitly of the Hamiltonian (19) 
with the linearized energy dispersion curves (23). From Eq. (33) we see 
that the Bose creation operators may be written as 



With the help of the Fermi anti-commutation relations we then find 

HQ,h'^^_ = hvF\q\b\, (87) 

where we have also made use of £p^k+q — £p,k = phvpq which holds for 
the linear spectrum (23). Now, a A^-particle ground state |A^4., A^_)o is 
an eigenstate of Hq with the energy 

Since |0, 0)o has zero energy by definition and the single particle energies 
(23) are counted relative to the Fermi energy, the result (88) is obtained 
readily by adding for Np > the single particle energies of the states 
occupied additionally in the iV-particle ground state |A''_|_, -/V_)o, while 
for Np < we have to subtract the (negative) energies of the particles 
removed. Using the commutator (87) we obtain for an arbitrary basis 
vector (39) 

(bir ( \ 

= \[^lL^\Ho + Y.nvF\q\mA\N+,N^)Q, (89) 

which shows that |A^+, A^-, {m-g}) is an eigenstate of Hq with eigenvalue 
EQ{N+,N^,{mg])=EQ{N+,N^)+hvFY, \q\mq. (90) 

^Soine authors use non-standard definitions of tiie Fermi annifiilation operator fields (45). 
To compare witli our notation, one has to make proper replacements, e.g., replace x by —x. 
In addition, the fields 'd(x) and (j>{x) are sometimes defined the other way round. 
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Since the basis (39) is complete and the basis vectors are clearly eigen- 
states of 



Hq = hvp 



97^0 P 



(91) 



with the proper eigenvalues (90), we see that Eq. (91) gives indeed a 
Bose representation of the Hamiltonian (19). 

Combining Eqs. (33) and (50), we may write the gradient of the 
Bosonic phase field ^p{x) as 



|:^-(") 



p 



V2^ 



E ^(m) 



i e*«^ ba + h.c. 



(92) 



q^O 



which gives 



L/2 



dx 



7^j:\i\M+bib, 



27r 



-L/2 



q^O 



(93) 



The Hamiltonian (91) may thus be written as 

L/2 



Hq = ttTivf I dx 22i 

-L/2 P L 



fd^ 



V dx 



^^ : +-liVp(iVp + l) 



(94) 



where : : puts the Bose operators in normal order. Further, in terms of 
the fields (58) this reads 



Ha 



flVp 



-L/2 '- P 



• (95) 



Since 11^ = —fid(f)/dx is the conjugate density to the phase field t? with 
the canonical commutator (62), we finally obtain the T?-representation 
of the Hamiltonian 



L/2 

Ho = ^ / dx 

-L/2 



1 



d^\ 



27r 



^^n|+U- :+72E^p(^P + l 



dx J 



L2 



(96) 



Likewise, using H^ = —hd'd/dx, we can readily write down a (^-re- 
presentation of Hq. 
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2.7. ACTION FUNCTIONAL 

In the usual way, we may introduce a Lagrangian 

L/2 



21 



Ln 



dxli,, 



d-d 



Hn 



(97) 



-L/2 



where the time rate of change of the T?-field reads 



l''4i^»-''i- 



OX n 



(98) 



which follows from the commutation relations (60). In the limit L ^ oo 
this gives 



Lo-- 


= h dx 


1 

2VF 


\dt) 


vf 
2 


\dx) 



(pI + /-) 



Nr, 



Pv 



(99) 



(100) 



is the average density of p-movers. Note that by definition the densities 
/Jp vanish in the ground state |0, 0)o, where the Fermi levels are at zt/cj?. 
The Pp determine a shift of the Fermi points, while the density fluctua- 
tions Pp{x) —'Pp arising from electron-hole pair excitations are described 
by the phase field ??. However, we may define a shifted phase field 



!?' = -(?+ ^/^T (p_|_ + p_) X - ^/vr (p_, 
with the properties 



p_) VFt 



(101) 



1 d^' 

/tF dx 

1 d-d' 



P+ + P- 



-VF (p^ 



(102) 



P- 



Now, the gradient determines the total density of right- and left -movers 
including the ground state density 'p_^_ +'p_ , while the time rate of change 
is proportional to the particle current. The Lagrangian then takes the 
simple form 



Ln 



Ijdx 


' 1 


/d^Y 


fd^Y 


VF 


[dtj 


^^ \dxj _ 



(103) 
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where we have omitted the prime. The noninteracting TL model with 
dispersionless spectrum (23) can thus be characterized by the classical 
action functional 



So = - I dt I dx 



±(«)^ .,(«)■ 

vp \ot J \dx J 



(104) 



which is the action of a harmonic string with wave velocity vp and a 
dimensionless displacement field d{x,t) measured in units of {Ti/vFpiY''^ 
where ^ is the mass density in the string. From this mechanical analogue 
it is obvious that Bosonization provides an alternative description of one- 
dimensional Fermions in terms of charge density oscillations rather than 
electron-hole pair excitations. 

We remark that in the general case of an arbitrary dispersion law 
Ep^k the Bosonization identities, in particular the representation (85) of 
the Fermi operators remain valid, however, the Hamiltonian Hq is no 
longer quadratic in the Bose operators bq, b]^. As a consequence, the 
mechanical analogue will be an anharmonic string which may, of course, 
be treated in the harmonic approximation when we restrict ourselves to 
low energy excitations. In combination with the Feynman path integral 
representation the classical action functional (104) can be a convenient 
starting point for quantum mechanical calculations. 

2.8. ELECTRON DENSITY OPERATOR 

When we use the TL model to make predictions for one-dimensional 
fermions, we have to remember that the model provides a local ap- 
proximation to the physical model in the vicinity of the Fermi points, 
while there is only a single energy dispersion curve for real electrons, cf. 
Fig. 1.1. This has consequences for observables like the density operator 

P{^) = 7 E e-('=-'=')- (4a,, - (4a,,)o) , (105) 

k,k' 

where the a,, a^, are Fermi operators for the underlying physical model 
with a single branch, and we have subtracted the constant density of the 
ground state with the Fermi points at zizkp. When we restrict ourselves 
to low energy states, the operators al.af^, can give a nonvanishing contri- 
bution only if both wave vectors k and k' are in the vicinity of the two 
Fermi points zizkp. However, for k ~ pkp, k' ^ p'kp, {p, p' = ±), the 
operators a]., a,, may be replaced by the Fermi operators cl,,, c ,,, of 
the TL model, and Eq. (105) splits into four terms 

pp' kk' 
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= ^(*t(^)*p,(a^) - (^iW^p'W)o) , (106) 

pp' 

where we have used Eq. (45) to obtain the second hne. Now, the diagonal 
terms {p = p') just give the densities (48) of p-movers, while for the non- 
diagonal terms we employ the representation (85) to find 

2iTa ^ ^ 

(107) 

where we have used the commutators (44) and the fact that ■d[x) and 

(j){x) commute at the same position x, which is seen from Eq. (60). 

Note that in the presence of a cutoff the sign function in Eq. (60) takes 

the form of the right hand side of Eq. (56) and therefore vanishes for 

vanishing argument. The density (106) may thus be written as 

p{x) = p+ (x) + p- (x) + p2kp (x) , (108) 

where 

P2. (x) = — (e'Hkp^+V^m+i (N.+N-.+l)] u\tj_ + h.C.) . (109) 

27ra V ^ / 

Hence, the density operator for real electrons is not just the sum of the 
densities in the two branches of the TL model, but there is an additional 
2fci;'-component p2kp {x) which comes from the fact that right- and left- 
movers are propagating in the same channel and interfere. 

In the limit L ^ oo with constant average densities (100), we may 
introduce the shifted phase field (101) with t = since the expression 
(108) gives the operator in the Schrodinger picture. We then have in 
view of Eq. (102) 

1 d 

p{x) = ^—d{x)+p2kA^), (110) 

vvr ox 



where we have again suppressed the prime on t9 which now contains the 
terms in the exponent of Eq. (109) that depend on the particle numbers. 
The operators U^U_„ in Eq. (109) associated with the scattering of an 
electron from the {—p)- into the p-branch can often be suppressed since 
changes of the particle numbers Np by 1 can be neglected for L ^ oo. 
Then p2kp (x) takes the simple form 

P2kp{x) = — cos [2kFX + 2^/^^{x)] , (111) 

TT 

where we have chosen a cutoff length a = kp of order a typical micro- 
scopic length. 
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2.9. 



FERMIONS WITH SPIN 



We now briefly summarize the modifications necessary to include the 
electron spin. Then, the electron spectrum has two branches s =T; i for 
spin up and spin down particles. For each species s, we may proceed 
exactly as for spinless fermions and define Bosonic phase fields t?s (x) and 
(ps {x) . There are now four branches of the TL model and correspondingly 
four particle number operators Np^s- When defining the ladder opera- 
tors, we have to include appropriate sign factors in the generalization of 
Eq. (40) to assure anti-commutation relations 



l'~^p,s ) '-^p', 



2 Sp^p> 6s^s' 



{U, 



p,s) 






2 5p,p' 5s,s' ( t^p,s 



(112) 



TJ Tr 



2 5r,.r,' (5, 



'p,p' ^s,s' 



that extend the relations (44) to the case of four branches. The Bose 
representation (85) of the Fermi operators takes again the same form for 
each species s, i.e., 



^pA^) 






It is often convenient to transform to the phase fields 



and 



^p(x) 

4)p{x) 



1 



— [d^{x)+^^{x)] 



^[0|(x)+</.^(x)] 



(113) 



(114) 



^a{x) 



(t>a{x) 



V2 



[M^)-^l{x)] 



-j=[(i)^{x) - 4>i{x)] 



that satisfy the commutation relations 

[^p(x),^p(x')]_ = K(x),^„(x')]. 



(115) 



Transport in Single Channel Quantum Wires 

[(l)p{x) , 0p(x')]_ = [0a (x) , 0ct(x')]_ = 



[^p{x),^„{x')]_ 



[0p(x), (t)a{x)\_ =0 
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(116) 



and 



[t9p(x),</<4x')]. 



[i?p(x) , 0p(x')] = [i?a(a:) , (l)a{x)\ = - sign(a; - x) . 



(117) 



These fields describe charge and spin density excitations. In particular, 
instead of Eq. (63) we now have 



2 d 

-^^p(a^) = P+a{x) + P+,i{x) + P^^{x) + P~,i{x) 

iV I t + iV , I + iV_ t + iV_ I 



"+,? 



"+,[ 



L 



(118) 



which gives the fluctuations of the total particle density of right- and 
left-movers, while 



2 d 

-^^a(a;) = p+,\{x) + /9_,|(x) - p+^i{x) - /?_4(x) 



7V+,I + A^_ 



TV. 



+,i 



iV_ 



L 



(119) 



determines fluctuations of the spin density. In the limit L ^ oo with 
given average densities 

(120) 



iV, 



p,s 



PiS 



L 



we may again introduce shifted phase flelds 'd'^ for each species s accord- 
ing to Eq. (101) and then obtain the classical action functional of the 
noninteracting model with dispersionless energy spectrum as 



So = ^T.J^tJdx 



1 



d^a\^ fd^r 



VF \ dt 



VF 



dx 



(121) 



where the sum is over the two spin directions a =t,i. Then, ^o is 
the sum of the actions (104) for each species s. The expression (121) 
remains, however, also valid if the sum is over a = p,a with the phase 
flelds introduced in Eqs. (114) and (115). Hence, the action functional 
for spinful electrons can be split in a charge and spin contribution 



'S'o — Sq p + Sq(j 



(122) 
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where each term has the form (104) of the action of a harmonic string. 
The representation of the true electron density operator (105) in terms 
of the Bosonic phase fields has also a straightforward extension to the 
spinful case. Here we give explicitly only the generalization of the ex- 
pression (110) valid in the limit L -^ oo. One now obtains 

(123) 





p{x) = 


\ TT OX 


02kF{x), 


where 




P2kF{x) = 


2kF 

= COS 

TT 


2kFX + \/27n?p(x) 


COS 


V27n?^(x) 



(124) 

is the 2/ci?-contribution from the interference between right- and left- 
movers. Again we have suppressed the ladder operators that need to be 
taken into account in general. 

3. INTERACTION, VOLTAGE BIAS, AND 

IMPURITIES IN THE 
TOMONAGA-LUTTINGER MODEL 

In the previous section we have investigated noninteracting Fermions 
in a single channel quantum wire and seen that the low energy physics 
of the Fermi gas can be described either in terms of occupation num- 
bers of single electron states or as excitations of Bosonic density waves. 
In higher dimensions noninteracting quasiparticles, supplemented by an 
electroneutrality constraint, give a rather accurate description also of the 
system in presence of Coulomb interaction, provided some parameters 
are replaced by effective parameters [24]. It has been known since quite 
some time that in one dimension the Fermi liquid description breaks 
down. Remarkably, the interaction leads to rather moderate modifica- 
tions of the action in the Bose representation. The real advantages of 
Bosonization will thus only emerge in this section where we first present 
the real, interacting TL model and then extend it to describe a quantum 
wire with a scatterer in presence of an applied voltage. 

3.1. ELECTRON-ELECTRON INTERACTION 

AND THE TOMONAGA-LUTTINGER 
MODEL 

We now take the interaction into account but restrict ourselves to 
spinless Fermions first. The electronic charge density in a quantum wire 
is largely compensated by a homogeneous positive background charge 
density. Therefore, we now fix the Fermi wave number kp so that the 
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state |A^+ = 0, N^ = 0)o, where all single particle states of p-movers 
with pk < kp are occupied, is electrically neutral. The density p{x) 
introduced in Eq. (105) multiplied by the electron change e is then the 
excess charge density in the wire, and the interaction may be written 

^int = \ j dx j dy p{x) U{x - y) p{y) . (125) 



The electron-electron interaction potential 



4 



U{x) = ^ E ^^ ^''^ (126) 

has real Fourier components Uq since U{—x) = U{x). Usually, the po- 
tential U{x) deviates from a simple Coulomb potential both at small 
and large distances. At small distances of order the lateral dimensions 
of the quantum wire, the wave function for the transversal motion of the 
strictly speaking three-dimensional particles becomes relevant. When 
these transversal components are integrated out, the resulting effective 
potential U{x) in the one-dimensional model remains finite for x ^ 0. 
On the other hand, at large distances one has to take into account the 
effect of gate electrodes and other nearby conductors that screen the 
long-range part of the Coulomb interaction. The effective potential U{x) 
then has a finite range i?, which implies that the logarithmic increase 
for small q of the Fourier transformed Coulomb potential is cut off at 
^-values of order R~^. In particular, C/g=o then remains finite, with the 
precise value depending on the geometry of the problem. We assume 
that this externally screened potential is still sufficiently long ranged so 
that Uq ^ U2kp- Then the 2/ci?-component of the electronic density 
(110) will give a negligible contribution when inserted in the interaction 
(125), and we find for the Bosonized interaction energy 

When we restrict ourselves to low energy excitations with wavelengths 
large compared to the range R of the interaction potential, we may re- 
place U{x) by a local interaction Uq 6{x), and then obtain the interaction 
term of the TL model 

*In the limit L ^ oo the sum ^ X/ '^ replaced hy ^ J dg . 
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We remark that a microscopic local interaction ~ 6{x) would of course 
have no effect on spinless Fermions as a consequence of the anti-commu- 
tation rules. Essentially, in Eq. (128) one neglects the wave number 
dependence of Uq for small q, but still U2kp '^ Uq. The interaction term 
(128) can readily be put in the action functional (104), which for the 
interacting model takes the form 



On 



h 



dt I dx 



1 

VF 



'dtj 



VF 



ttUvf) \dxj 



(129) 



and thus remains of the form of the action of a harmonic string. This 
clearly shows the great advantage of the Bose representation: We still 
have a model of free Bosonic charge density excitations, only the wave 
velocity is altered by the interaction, while in the Fermi representation 
the interaction leads to quartic terms in the Fermi operators. 
It is customary to introduce the coupling constant 



9 



1 + 



^0 \~'/' 



irhvF 
and the charge density wave velocity 

- }!L 
9 



(130) 



fl31) 



In terms of these quantities the action functional of the TL model reads 



Sp = — I dt I dx 



'dt 



dx J 



(132) 



It would go beyond the scope of this article to demonstrate that the 
action (132) indeed describes the low energy properties of spinless one- 
dimensional Fermions correctly. Here we refer to the literature [8, 15]. 
We would, however, like to point out that the matter is in fact more 
complex than what our plausible "derivation" of Eq. (132) might sug- 
gest. The Coulomb interaction is strong and affects all states not only 
those near the two Fermi points. Hence, it needs to be introduced in 
the underlying physical model with a single energy dispersion curve £k- 
Afterwards, one may integrate out states far from the Fermi points until 
one reaches energy scales sufficiently close to the Fermi energy to allow 
for a linearization of the spectrum. At this point the Hamiltonian may 
be re-written in terms of the right- and left -movers of the TL model, 
but the parameters of the model are then already renormalized by the 
aforementioned elimination of high energy excitations and additional in- 
teraction vertices are generated. One can, however, conclude from a 
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renormalization group study [7] that the action (132) is indeed a low en- 
ergy fixpoint of a (spin-polarized) quantum wire. From these remarks it 
is clear that the parameter C/q ™- Eq. (128) does not necessarily coincide 
with the Fourier coefficient at g = of the interaction potential. This 
latter quantity can be considered an estimate of Uq which becomes more 
accurate for large electron densities at which the effect of the Coulomb 
interaction is weaker. In the sequel we use g and v as fundamental pa- 
rameters of the model. For repulsive interaction C/q > and thus g <1. 

3.2. SCREENING OF EXTERNAL CHARGES 

Let us assume that we perturb the quantum wire by an external charge 
density 

/3cxt(a:,t) = eQ{x,t) , (133) 

which interacts with the electronic charge density ep{x, t) via the same 
effective potential U{x) introduced in the previous section. We are inter- 
ested in the long wavelength response of the quantum wire and may thus 
disregard the 2/ci;'-component of the electronic charge density. With the 
local approximation Uo5{x), the action (132) is then modified to read 

S = Sp - ^ f dt f dxQ(x,t)-^^(x,t). (134) 

^J^T J J OX 



Since the action has only terms linear and quadratic in t?, the average 
electronic density {p{x, t)) caused by Q{x, t) can be determined from the 
phase field -dlxji) minimizing the action (134). The equation of motion 

obeyed by the minimal action field is readily solved in terms of the 
Fourier representation 

^(x, t) = ^ Jdqjduj ^q, u) e*-?-^-* . (136) 



We find 



5/ N gUo ivqQ{q,uj) 

^q,uj) = ^-^ 2"^, (137) 



which yields for the electronic density {p{x, t)) = -j= ^ 'd{x, t) in Fourier 
space 

{p[Q,^)) = [^- 9 )^ 2"!-, (138 
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where we have expressed Uq in terms of g and v by means of Eqs. (130) 
and (131). 

Now, the relation between the external charge density eQ{x, t) and the 
resulting screening charge density e{p{x,t)) is governed by the dielectric 
function 

e{g,u.) = ^^^^^ . (139) 

Q{q,uj) + {p{q,uj)) 

Combining this with Eq. (138) we find 

uP' — iP'cP' 
e{q,uj) = — ^-^. (140) 

In particular, in the static case w = we have 

e = \ = l + ^, (141) 

which shows that the interaction parameter g is directly related to the 
dielectric constant of the quantum wire. In a metallic system the dielec- 
tric function has a pole for to ^ 0, q ^ associated with the perfect 
screening of static charges leading to electroneutrality. However, in the 
TL model the long range part of the Coulomb interaction is assumed to 
be screened by other conductors as explained above. Then Uq is finite 
and there is a finite dielectric constant in the zero frequency and long 
wavelength limits. 

This is in accordance with the fact that the total screening charge in 
units of e 

Q,= fdx{p{x)) = {piq = 0)) (142) 



accumulated near a static impurity charge Q follows from Eq. (138) as 
[25] 

Qs = -{l-9^)Q- (143) 

Hence a fraction g^Q of the external charge remains unscreened, and 
the quantum wire is in general not electroneutral. As we will discuss 
in greater detail in the next section, the charge g^Q is screened by the 
electrode responsible for the finite range of the interaction. Formally, 
the limit of long range Coulomb interaction corresponds to g ^ which 
implies electroneutrality of the wire. 

We mention that apart from the long wavelength response of the quan- 
tum wire to an impurity charge there is also a 2/ci7-response leading to 
Friedel oscillations of the charge density. We will not discuss this here 
but refer to the recent literature [26, 27]. 
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Figure 1.6. Electrostatic model of a quantum wire coupled by a capacitance per 
unit length cq to a gate. In a nonequilibrium wire the band bottom is shifted by an 
electric potential tp. 



3.3. ELECTROSTATICS OF A QUANTUM 

WIRE 

As we have seen in the preceeding section, the electrode screening the 
long range part of the Coulomb interaction plays an important role in 
the electrostatic response of a quantum wire to external charges. We can 
visualize the TL model as a one-dimensional quantum wire screened by 
a gate coupled to the wire by a distributed capacitance as depicted in 
Fig. 1.6. The interaction energy (128) can then be interpreted as the 
charging energy of the wire-gate capacitance 



H, 



Uo 



int 



dx p{x) 



dx 



e'^p{x) 



2 j "■ ' J 2co ' 

where the capacitance per unit length cq is determined by 



Co 



ttUvf ( -5- 
9^ 



1 



(144) 



(145) 



An electronic charge density ep{x) in the wire polarizes the capacitance, 
and the resulting electric potential 



ip{x) 



ep{x) 
Co 



(146) 



shifts the band bottom of the quantum wire. This is directly related 
to the underscreening of external charges discussed previously as can 
be seen from the following consideration. In a noninteracting system 
an increase of the Fermi energy by AE shifts the wave number of the 
Fermi points from zizkp to i{kp + Ak) where Afc = AE/fivp. There 
are (L/27r)Afc single particle states in the wave number interval Ak and 
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thus the density of p-uioveis increases by 

Here the factor 1/2tt%vf gives the density of states of p-movers at the 
Fermi energy. Accordingly, the electronic density increases by 

A P 

A/9 = A/9+ + Ap_ = — — . (148) 

ttTivf 

On the other hand, in an interacting, electroneutral system the shift of 
the Fermi energy by /S.E is accompanied by a shift of the band bottom 
by the same amount and the electronic density remains unchanged. In 
the TL model the situation is in between these two extremal cases. The 
band bottom is shifted by eip where (/? is the electric potential differ- 
ence between the wire and the gate electrode. Then, the change of the 
electronic density is determined by 

AP=^^. (149) 

In view of the relations (145) and (146) we have 

e2 / 1 \ 

ev? = — A/9 = TT^fF ( — - 1) A/9, (150) 



Co \9 



which gives 



Ap=^——. 151 

ttTivf 

Hence, we see again that a fraction g^ of the "bare" charge density (148) 
caused by the shift of the Fermi level persists as true electronic charge 
density in the wire. 

3.4. VOLTAGE BIAS AND BOUNDARY 

CONDITIONS 

It is not difficult to generalize the preceeding considerations to a 
nonequilibrium quantum wire in presence of an applied voltage. Let 
us consider a quantum wire which is attached at the ends to two- or 
three-dimensional Fermi liquid reservoirs. We assume that the contacts 
between the wire and the reservoirs are adiabatic, which means that 
at the ends the quantum wire widens sufficiently slowly to avoid any 
backscattering of outgoing particles into the wire. This is the usual 
assumption underlying Landauer's approach [16] to the conductance of 
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mesoscopic wires. If the contacts are not adiabatic, there will be an ad- 
ditional resistance depending on the precise realization of the contacts 
and not only on intrinsic properties of the quantum wire. 

In equilibrium there is an equal amount of right- and left-movers in 
the wire and the electrochemical potential is constant. When we attach 
the wire to reservoirs, the influx of right-movers at the left end of the 
wire will depend on the electrochemical potential of the left electrode 
which we assume to be eUi above the Fermi energy of the equilibrium 
quantum wire. Ul is then the voltage between the left reservoir and the 
gate electrode screening the wire. In the absence of interactions, the shift 
of the Fermi energy by cUl would increase the density of right -movers 
near the left end of the wire by 

bare _ ^^L (152) 

2-KhvF 

as depicted schematically in Fig. (1.7) a. Note that below we will con- 
sider quantum wires with impurities. Then, the reservoir determines the 
density of incoming particles only in the clean section of the wire near 
the end, where the incoming particles have not yet interacted with the 
impurities. The density of outflowing particles, on the other hand, will 
be affected by impurities. 

In the presence of Coulomb interaction the excess charge density 
caused by the reservoirs will charge the distributed wire-gate capaci- 
tance leading to a shift of the band bottom by 

e(f = — p = nhvF -;r — 1 p , (153) 

Co Vr / 

where ep is the true charge density in the wire that has to be determined 
selfconsistently. The true density of right-movers near the left end of 
the wire is then, cf. Fig. (1.7)b, 

e{UL-V>) ..... 

which means that the bare charge density (152) is partially screened. 

In terms of the Bosonic phase field t?, we have by virtue of the relations 
(102) 

1 fdd 1 d{)\ ,,,,, 

where we have omitted the prime on •& as in all equations following (102). 
We remark that in the relation 

1 (9i9 
-^— = -vf{p+ - p~) (156) 
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Figure 1.7. Density of right -movers in a quantum wire attached at the left end 
to a reservoir with voltage Ul against the gate electrode. The equilibrium density is 
represented as light-grey area and the nonequilibrium excess density as dark-grey area, 
a) Sketch of the situation in the absence of Coulomb effects and b) with interaction 
taking the shift etp of the band bottom into account. 




Transport in Single Channel Quantum Wires 35 

the Fermi velocity vp is not altered by the Coulomb interaction such as 
the velocity of charge density waves, because the microscopic expression 
of the particle current operator in terms of Fermi operators is indepen- 
dent of the interaction in the absence of a vector potential [28]. For the 
electric potential (p we obtain from Eqs. (153) and (102) 

/I \ d'd 

e(^ = ^Tit;^ f — - Ij — , (157) 

which can be combined with the relations (154) and (155) to give for x 
near the left end of the wire 

(158) 



g'^ dx vp dt ^^hvp 

Similar considerations hold for the density of left-movers near the right 
end of the wire attached to an electrode at potential Ur. Strictly speak- 
ing, all of the above considerations hold for the average densities of right- 
and left-movers. Hence, we find that the coupling to reservoirs can be 
described in the Bosonized TL model in terms of radiative boundary 
conditions for the phase field [17] 



g"^ dx Vp dt) ' ^ 2 ^/^Thvp ' 



Id ^ d\ .^. .. eUn 



(159) 






y/T:hvp 



If we impose these conditions at a point near the reservoirs, they will 
be seen to be obeyed at any point in the impurity-free clean sections at 
the ends of the wire. 

3.5. IMPURITY POTENTIAL 

Nontrivial dc transport properties of the quantum wire do not arise 
from the Coulomb interaction and the resulting non-Fermi liquid behav- 
ior alone but only in connection with impurities causing backscattering of 
electrons. A nonmagnetic impurity at position xq couples to the charge 
density and gives rise to an energy 

Hw = I dx W{x - xo) p{x) , (160) 

where W{x) is the impurity potential. When we want to add this cou- 
pling term to the TL model, we have to note that the perturbation af- 
fects all states, also those far from the Fermi points, and we should again 
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start from the underlying physical model with a single branch. Applying 
essentially the same line of reasoning used to include electron-electron 
interactions in the TL model, we obtain for the case of an interaction 
with an impurity at position xq 

Wo d^ 



F, 



w 



IT dx 



+ 



^^^^^^ cos [2kFXo + 2^i?(xo)] . (161) 



Of course, this form results when we insert the representation (110) of the 
electronic density into Eq. (160) and then replace the Fourier coefficients 
Wq of the impurity potential by constants for g ~ and \q\ ~ ^kp, 
respectively. Since an elimination of states far from the Fermi points is 
necessary before we can rewrite the Hamiltonian in terms of the Fermi 
operators c, 



p,A;' p,k 



of the TL model, the coefficients Wq and W2kp must 
be interpreted as effective parameters that are not directly related to 
Fourier coefficients of the microscopic interaction potential. In addition, 
the elimination of states will generate higher order processes involving 
several electrons with momentum transfer ikp, Gkp, and so on. It turns 
out that the 2A;i7-processes dominate at low temperatures [18]. 

Although an impurity is typically laterally displaced from the center of 
the wire, the potential W{x) will essentially have the same behavior as 
the electron-electron potential U{x) discussed previously. Thus, from 
the bare Fourier components we would conclude W2kF *^ Wq. But, 
as we shall see, the forward scattering Wq, which scatters a p-mover 
into another state of the same branch, has no effect on the transport 
properties of the wire, while even a small backscattering term W2kp has 
a dramatic effect at low energy scales [18]. 

To show that forward scattering is unimportant, we write the TL 
Hamiltonian in presence of an impurity in the form 



H 



hvp 



dx 



dx 



+ 



dx J 



+ 



2tt 



dx 



dxj 



(162) 



+ 



Wa d^ 



dx 



+ ^^^^hE^lL cos [2kFXo + 2V^T?(xo) 



vr 



which includes the Hamiltonian (95) of the noninteracting model, the 
interaction (128), and the impurity term (161). Now, the unitary trans- 
formation 



U = exp 



shifts the charge density, since 



i-v/vr / dxa{x) 4>{x) 



U^-U = T- + Vvra(x) 
dx dx 



(163) 



(164) 
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This is readily shown by considering the auxiliary function 



Fis) 



isA ^_^ „-isA 

ox 
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(165) 



where A = y/ir ^ dx a{x)(l){x) is the exponent oiU. One then finds 
d 



ds 



F{s) = e'''^i[A, d{)/dx\_ e 



isA 



TT a[x] 



(166) 



where the commutator is evaluated by means of Eq. (60). Since -F(O) = 
d-d/dx, we find -F(l) = d-d / dx + ^/tt a{x) which is just the relation (164). 
Based on the same commutator, one also finds 



where 



7/ 2J0ri?(a;)7/-l _ 2i^/ni}{x)+iri(x) 



rj{x) = IT dy sign(x — y) a{y) 



(167) 



(168) 



To see this, we write the left hand side of Eq. (167) as e* e* e~ with 
the operator A introduced previously and B = 2-^/7rT?(x). Then, using 
twice the relation (83), we find 



e'^ e'^ e~'^ 



^iA ^i(B-A) g-i[A,B]_ ^ giB-[A,B] 



(169) 



which gives the transformation (167). 

With the help of the relations (164) and (167), the Hamiltonian (162) 
is transformed into 






ox J g'^ \0x J 



\pK dx 



+ 



X=Xt) 



9' 



/ dxaix) — — 
J ox 



(170) 



H — cos {ikpXQ + 2^/7rT?(xo) + ??(xo)] 



vr 



where we have omitted terms that depend only on a{x) but not on 
the phase fields. Further, we have introduced the interaction parameter 
(130). Now, with the choice 



q(x) 



Trhvp 



S{x - xo) 



(171) 
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the terms in the second Hne of Eq. (170) cancel, and we are left with the 
Hamiltonian 



H' 



flVp 



+ 



2 

W2kpkF 



dx 



'W 



dx J g"^ \dx J 



d^Y 



(172) 



TT 



COS [2kFXo + 2^/TT'd{xo) 



which contains only a backscattering term. Note that for an impurity 
with charge eQ giving rise to the Coulomb potential W{x) = QU{x), 
the quantity 



a(x) 



Trhvp 



Q 6{x - xo) = -(1 - g'^)Q S{x - xq) 



(173) 



is just the screening charge density (138) caused by the static impurity 
charge. The unitary transformation thus removes the screening cloud, 
which is the main effect of the forward scattering term. Below we will 
determine the current-voltage relation of a quantum wire in presence of 
a single impurity at position x = 0. This study will be based on the 
action 



s = A 

2s 



dt / dx 



1 fd^^^ 
V \~dt 



'df\' 
dx) 



X I dt cos [2^/tt'&{x = 0, t)] 



(174) 



associated with the Hamiltonian (172), where A characterizes the impu- 
rity strength. 

3.6. INTERACTING ELECTRONS WITH 

SPIN 

Here we briefly summarize the changes necessary to include the elec- 
tron spin. For noninteracting electrons we found that the action of spin- 
ful electrons can be split into a charge and a spin contribution. Since 
the Coulomb interaction couples only to the charge density, we might 
argue that the charge part of the action is modified by the interaction 
in the same way as for spinless electrons while the spin part remains 
unchanged. However, this would ignore the fact that the Coulomb in- 
teraction must be introduced in the underlying physical model and the 
transcription to TL Fermions can only be made after an elimination of 
high energy excitations. This has consequences also for the spin density 
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waves, in particular, the spin wave velocity Vfj becomes smaller than the 
Fermi velocity vp ([29]). The low energy physics is then governed by the 
action 



s = i- 

2s 



dt i dx 



dt I dx 



v\ dt 



1 fddr. 



dt 



dx 



dx 



(175) 



The TL model is thus characterized by three parameters v, v^ and g. 
There is no coupling constant §„ for the spin sector, which can be traced 
back to spin rotation invariance [8, 15]. A detailed discussion of these 
issues would go beyond the scope of this article and we refer to the 
literature cited. As in the noninteracting case, the action splits into a 
charge and a spin part. However, the difference between the charge and 
spin wave velocities in an interacting wire has important consequences 
and leads to the notable phenomenon of spin-charge separation. When 
an electron is injected into a wire it causes a charge and a spin pulse 
propagating with different velocities. 

In case an impurity is added at position x = to the model we see 
from the 2A:F-part of the electron density in the spinful case (124) that 
the important backscattering term now has the form 



X dt cos 



2TT^p{x 



0,t) 



COS 



2TT^„{x = 0,t) 



(176) 



Hence, the impurity couples the charge and spin sectors making the the- 
ory of dirty quantum wires in the absence of a spin-polarizing magnetic 
field more involved. 



4. CURRENT- VOLTAGE RELATION OF A 

QUANTUM WIRE 

In this section we apply the theory developed so far and determine 
the current in a quantum wire with an impurity as a function of the 
applied voltage and the temperature. Rather than giving an overview 
of the results available in the literature, we treat a special case in some 
detail to illustrate how the formalism explained in the previous sections 
can be employed to obtain concrete results. With this background read- 
ers should then be prepared to embark on reading the recent original 
literature on the subject. 
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Figure 1.8. Gated quantum wire adiabatically connected to reservoirs at voltages 
Ul and Ur. a backscatterer is located in the center of the wire. 

4.1. PARTICULAR SOLUTION AND 

FOUR-TERMINAL VOLTAGE 

We study a single channel quantum wire with an impurity at position 
X = 0. The conductor is attached to reservoirs at voltages Ul and Ur 
relative to the gate electrode screening the wire as sketched in Fig. (1.8). 
The applied voltage 

U = Ul-Ur (177) 

will then drive a current / through the wire. For simplicity, we shall 
restrict ourselves to the case of spinless electrons. We can then base the 
consideration on the action functional (174) and the boundary conditions 
(159). 

As discussed previously, the phase field '&{x^ t) describes also the mean 
particle densities (100) due to a shift introduced in Eq. (101). To take 
care of the boundary conditions, we look for a particular field 6{x,t) 
that satisfies the wave equation in the clean parts of the wire and the 
boundary conditions. This solution is of the form of the shift in Eq. (101) 

2 

e{x,t)= ^l [(Ul + Ur)x-V\x\]- -^{UL-UR-V)t, (178) 
2^/Trhvp Zy/TTh 

where V is an arbitrary parameter. We now split the phase field '!?(x, t) 
into 

'd{x,t) = e{x,t)+ip{x,t). (179) 

The deviation ip{x,t) from the particular solution (178) will then obey 
equilibrium boundary conditions. Further, we fix the parameter V by 
the requirement that in the steady nonequilibrium state 

^((^(x,i))=0. (180) 
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Now, the average current I may be written as 

/ = et;^(pH_-p_) = --^ -(.?), (181) 



where we have used Eq. (102) to obtain the second equahty. Hence, the 
condition (180) means that the average current is determined solely by 
the particular solution (178), and we have 

I=^-{U-V). (182) 

Since in a steady state the current / is independent of x, we may impose 
the condition (180) fixing the voltage V at any point x. 

The particular solution (178) not only determines the average current 
but also the average charge density in the nonequilibrium quantum wire. 
Since {ip) obeys equilibrium boundary conditions as well as the condition 
(180), it does not contribute to the average density (110) 



where we have omitted the 2A;i?-component which gives an additional 
oscillatory contribution near the impurity. This Friedel oscillation com- 
ponent is not seen in a density smoothed over length scales of order 
\f = Itx jkp. From the particular solution (178) we obtain for the aver- 
age density 

The first term just describes the change (151) of the average electronic 
density as a consequence of the average shift e(C/i + Ur)/2 of the Fermi 
energy. This term is absent if the voltage U is applied asymmetrically, 
i.e., for Ul = —Ur = U/2. The second term gives an asymmetric 
component of the charge density in presence of an impurity. The density 
drop 

AP=4^ (185) 

across the impurity site is associated with a difference 

A^ = g^eV (186) 

of the effective chemical potential on both sides of the impurity. Fur- 
thermore, the drop of the charge density across the impurity site is also 
associated with a change of the electric potential (153) by 

A,. = ^(l-l)Ap=(l-,^)y. (187) 
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In a hypothetical ideal measurement of the voltage drop across the impu- 
rity site one would observe the difference of the electrochemical potential 
[16] 

Acp + -An = V. (188) 

e 

Hence, the parameter V introduced above coincides with the average 
four-terminal voltage V, which is the part of the applied voltage U drop- 
ping across the scatterer.^ 

In view of Eq. (182) the determination of the current-voltage relation 
corresponds to a calculation of the four-terminal voltage. Two limiting 
cases are evident from physical grounds. In the absence of a backscat- 
terer (A — > 0) we have V = and obtain from Eq. (182) 

I = GoU, (189) 

where Go = e^ /h is the conductance of a clean wire. This is the same 
result as obtained previously for noninteracting electrons. Hence, for a 
clean quantum wire with adiabatic contacts to the reservoirs the inter- 
action has no effect on the conductance [30, 31, 32, 17]. On the other 
hand, for a very strong backscatterer (A -^ oo) we have V = U and the 
current / vanishes. In the remainder we shall discuss how we get from 
one limit to the other. 

4.2. PATH INTEGRAL ON THE KELDYSH 

CONTOUR 

To treat the nonequilibrium quantum wire with the action functional 
(174) quantum mechanically, we have to evaluate a Feynman path inte- 
gral on the Keldysh contour. For an introduction to the Keldysh tech- 
nique we refer to the review article [33], however, the basic idea can be 
understood in the following way. Assume that at time to the system is 
described by the density matrix VF(to) and let H be the Hamiltonian 
including the coupling to the reservoirs. The density matrix at a later 
time tf is then given by 

W{tf) = e-i^(*/-*«)) W{to) eS^(*/-*°) . (190) 

Each of the two time evolution operators e^'i^'*-''^*''^ may be written as a 
Feynman path integral. Since we are interested in steady state properties 
independent of the initial state VF(to), we take the limit to ~^ — oo. The 

^ The readers should be aware of the fact that in some early treatments of transport properties 
of the TL model the discrimination between U and V was not made. 
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Figure 1.9. The Keldysh contour runs from — oo to tf along the real time axis 
and back to — oo. The two branches of the contour correspond to the forward and 
backward propagators in Eq. (190). 



trace over such a time propagated operator leads to a path integral of 
the form 

Zo= /"p[^]e^^[''l, (191) 

where S['d] is the action functional (174) with the time integration /dt 
running along the Keldysh contour depicted in Fig. (1.9). In order to 
employ this path integral for the calculation of expectation values, we 
first decompose the phase field according to Eq. (179). The action (174) 
then takes the form 



— / dt 
25 



dx 
dtip{0,t) 
X I dt cos 



dt) 



dip 
dx 



eV_ 
/i 



(192) 



2V^ip{0,t) --{U-V)t 



where we have omitted terms independent of ip{x.,t). Further, we have 
taken into account that 

dt^=0 (193) 



for a time integral along the Keldysh contour and 

d 



f d 

/ dx sign(x)— — (p{x, t) = — 2(/3(0, t) . 



(194) 



For the phase field (p{x, t), which obeys equilibrium boundary conditions, 
the nonequilibrium situation becomes apparent in two modifications of 
the action (192). The term in the second line comes from the voltage 
drop across the impurity site. It describes the potential energy 



dx 



-— sign(x) 



e dip{x,t) 
IT dx 



eV , , 
^(/.(0,i) 



(195) 



of a charge density fluctuation (e/y^) -^ in presence of an electrochemi- 
cal potential —iy 12) sign(x). The remainder U—V of the applied voltage 
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U shows up in the third hne of Eq. (192) as a Josephson-type phase shift 
in the pinning potential caused by the impurity. 

Instead of the path integral (191) we now study the functional 



Z[ri]= fv[^]e^^^^^+'^S'^^ 



r,{t) |-</p(0,t) 



(196) 



where S[f] is the action (192), and where we have introduced an aux- 
iliary field r]{t) on the Keldysh contour. Z[r/] is a generating functional 
for expectation values of §f:f{0,t). In particular, the condition (180) is 
now equivalent to 



5r,{t) 



0, 



»?=o 



(197) 



which needs to be evaluated to determine the four-terminal voltage V. 



4.3. 



COULOMB GAS REPRESENTATION 



In the sequel we outline one of the methods available to determine the 
current-voltage relation. While for the one impurity problem considered 
here, an approach based on the thermodynamic Bethe ansatz is most 
powerful [34], we present here a technique which remains useful also for 
multi impurity problems. First, we split the Keldysh contour explicitly 
into the two branches and denote the phase field on the branch from 
— cx) to tf by ip{x, t) and the field on the branch from tf back to — oo by 
ip'{x,t). Further, we split the action (192) into two terms 



s„ = A 



2gJ- 



dt / dx 



eV 



dt J \dx J V 



dt [(^(0,t)-(^'(0,t)] 



dt 



+ v 



dip' 
dx 



(198) 



and 



Sx = —X I dt { cos 
— cos 



2^p{Q,t)--{U -V)t 
n 

2^^'{0,t)-^{U-V)t 



With the help of the trigonometric relation 

cos a — cos (3 - 



a + (5 . a — (3 
2 sm sm 



(199) 



(200) 
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the action S\ may be written as 

Sx = 2\ r dt cos A{t) sin B{t) , (201) 

J —oo 

where 

A{t) = V^[^iO,t)+ip'iO,t)] -^{U-V)t + 5 

(202) 
B{t) = ^[ip{0,t)-ip'{0,t)]. 

From Eqs. (199) and (200) we obtain 5 = — |, however, the precise value 
of this phase must be irrelevant, since we can always add a constant 
phase to the particular solution (178), e.g., by replacing t by t — to- Such 
an additional phase of 9{x, t) leads to a shift of 5. 

Using Eq. (201) we find by expanding in powers of A 



e^"^^ 



1 + Y. ( X j J ^"* n cos A(tj) sinS(tj) , (203) 



where we have introduced the abbreviation 



Vnt = I dtn I dtn~l • • • / dti . (204) 



oo •! —oo J ~oo 



Next, we write the trigonometric functions as 
cos A{tj) = 7j J2 



2 ' 



(205) 



1 



sin B{t^) = - E 



2i 
which gives 

-rV At;i 



Vje 



ivjB(tj) 



"^ = 1 + E E [U^] J^nte^.^^, (206) 






where 



1 " 

-Sn = Y. i^jMtj) + VjB{tj)] (207) 

is linear in the phase fields ip, (p' by virtue of Eq. (202). The benefit of this 
expansion is that now the path integral (196) is Gaussian order by order. 
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Therefore, we can integrate out the (p and (p' fields. Essentially, the 
calculation of these path integrals amounts to an explicit computation 
of the fields minimizing the action. 

On a formal level we may consider the auxiliary variables Uj, Vj as 
charges on the time axis that are coupled to a harmonic string described 
by the phase fields. The integration over the fields c^, ip' then corresponds 
to the elimination of a harmonic bath in the context of dissipative quan- 
tum mechanics. This problem is well studied in the literature [35, 36, 37]. 
It would go beyond the scope of this article to present these methods 
explicitly here. Once the Gaussian integrals over the phase fields are 
carried out, the sum over the charges Vj can be done straightforwardly, 
and the generating functional (196) is obtained in the form 



Z[r]\ = e^^{- I dt I dt' 7]{t) C{t - t') 7]{t') - igj^V j dt 7]{t) 



(208) 
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dt 7]{t)C{t - tj) -i-{U-V + gV)tj 



2m- 1 

sin [ngr]{t2m)] H ^^^ ( ^^ 



2m 
k=j+l 



(209) 



Here the sine functions arise from the sum over the charges Vj by means 
ofEq. (205). The function 



C{t) = 2g\n 



/3A . (T,\t\ 



-K 



(210) 



describes an effective interaction between the charges Uj whereby /3 = 
l/ksT is the inverse temperature of the quantum wire and A is the 
cutoff energy of the harmonic string. The temperature emerges from the 
asymptotic conditions on the fluctuations of the phase field for t ^ — oo. 
To simplify notation we have chosen the same auxiliary field r]{t) on both 
branches of the Keldysh contour so that the second term in the exponent 
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of (196) reads ^/vr j'J^dtrjit) ^ [99(0, t) + v3'(0, t)]. Furthermore, we have 
taken the hmit tf ^ 00. 

In the series (208) only even terms in A survive, because C{t) is a 
long range interaction which suppresses all terms that do not satisfy the 
charge neutrality condition 



E 



Ui 



0. 



(211) 



Since the charges Uj = ±1, this condition can only hold for n = 2m. 
The constraint (211) is indicated as a prime at the sum over the Uj in 
Eq. (209). Because of this condition the phase 6 in Eq. (202) drops 
out. The representation (208) of the generating functional is known 
as the Coulomb gas representation of the problem, since some terms 
allude to the partition function of one-dimensional charges interacting 
with the "Coulomb potential" (210). Note that in view of the factors 
sin[7r(7?7(t2m)] the generating functional (208) obeys the normalization 
Z[r] = 0] = 0. 

The four-terminal voltage V can now be determined from the con- 
dition (197). A nonvanishing contribution of Zm[r]] only arises if the 
variational derivative acts upon sin[7r(777(t2m,)]. Introducing the time- 
difference variables 



we find 



where 



h 



V = K(-[U~V + gV] 



(213) 



Kin) = vrim ^ (-1)™ f-j / dn ... / dT2m-i (214) 

m=l 



X Y^ exp 



2m 1^^^ 

j>k=l \l=k 



2m- 1 

n sin(7r5(pj) 



-ipjflTj 



Here 



2m j 

Pj = E Uk = -^ Uk 

k=j+l fe=l 



(215) 



is the accumulated charge in the time interval t > tj. Sums of the 
form (214) with an interaction (210) are again familiar from dissipative 
quantum mechanics [38] 
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4.4. EXACT SOLUTION FOR G = | 

The explicit evaluation of the sum (214) is greatly simplified in the 
special case g = -^^ since the sine-functions then suppress many terms. 
The so-called "collapsed blip approximation" [38] of the Coulomb gas 
problem becomes exact for 9=2 ^^'^ K{Q) can easily be calculated 
explicitly. One finds 

where ip{z) is the digamma function and 

Xb = n^ (217) 

is a renormalized energy scale of the impurity. The four-terminal voltage 
V then follows from Eq. (213), which for g = 2 reads 

^V = Ki^^[U-V/2]^ . (218) 

It can be easily seen that the solution is of the form 

eV _ (eU kBT\ 
Xb \Xb ' Xb J 

This shows that the solution exhibits scaling. When the energies eU and 
ksT are measured in units of the renormalized impurity energy scale A^, 
we find the same four-terminal voltage (in units of A^) for any impurity 
and hence also the same form of the current -voltage relation (182). 

In Fig. (1.10) we depict the current-voltage curve for various tem- 
peratures. We see that the interaction has a dramatic effect on the 
/-^-characteristics in presence of a scatterer. While in the noninteract- 
ing model the conductance is essentially independent of /c^T and eU for 
small energy scales, there is a strong suppression of the conductance in 
the interacting system for energies below A^. At zero temperature 



I _ 1 (eU\ 
To ~ 24 [x^J 



2 
for eU « Xb (220) 



and 



— = 1 - TT^ for et/ > Ab , (221) 

Jo ell 

where /q = {e^ /h)U . Hence, the differential conductance shows a perfect 
zero bias anomaly and vanishes ~ C/^ for U ^ Q. On the other hand, the 
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Figure 1.10. Current-voltage characteristics for several temperatures T. The cur- 
rent is normalized to 7o = [e^ /h)U . 



conductance approaches the value e^ /h of a clean wire for large voltages. 
In practice, this latter limit will not be reached, since for large U the 
low energy description in terms of the TL model breaks down. On the 
other hand for T > the linear conductance for [/ ^ reads 



G{T) = - 



Q-^'{ 



+ h: 



h Q + ifj' 



:\+h) 



(222) 



where = 2T:kBT /\b- Hence, the conductance vanishes ~ T^ for 
ksT <^ Xb- We mention that the exact solution for g = ^ was first 
obtained by re-Fermionization techniques [17]. 



4.5. 



CONCLUDING REMARKS 



The explicit results for the cur rent- volt age characteristics derived in 
the preceeding section show that for a spin-polarized quantum wire the 
crossover from the weak impurity problem at higher temperatures or 
larger applied voltages to the strong coupling problem at low tempera- 
tures and small voltages can be solved exactly. The scaling function can 
also be obtained for arbitrary g < 1 [34]. Mostly, one has to be content 
with more limited results. In the spinful case the current is known only 
asymptotically for large and small energy scales [18]. Several topics are 
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still under investigation. For instance, in the presence of many impurities 
additional features such as resonant tunneling may arise. Another area 
of active research is the current noise in one-dimensional conductors. 

We have focussed here on the physics of screened single channel quan- 
tum wires at low energy scales where the finite range of the interaction 
is unimportant. Many results can in fact be extended to finite range 
interactions. However, in the absence of a gate, the long range nature of 
the Coulomb interaction requires special attention [39]. There are many 
similarities between the transport properties of quantum wires and edge 
currents in fractional quantum Hall bars [40]. Since in these latter de- 
vices right- and left -moving branches are spatially separated, they are 
described by a chiral TL model [41, 42, 43] with somewhat different 
transport properties. The methods presented here are also useful to 
model the low energy electronic properties of carbon nanotubes [44, 45]. 
Again, due to features of the bandstructure, some differences arise lead- 
ing to four channels where one channel is charged as in a semiconductor 
quantum wire. We hope that this article will facilitate the study of 
these and many other exciting new results on transport properties of 
one-dimensional Fermionic systems. 
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